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1.  INTRODUCTION 


In  the  computation  of  electron  residence  time  within  the  core  cavity  of 
Folywell^"'/HEPS  devices,  it  is  necessary  to  determine  the  average  time  of  transit  of  the 
electrons  in  the  system  across  the  system.  This  is  needed  in  the  EKXL  code  as  a  parameter 
in  the  basic  code  algorithm  rdating  electron  edge  densities  at  successive  time  increments 
within  each  time  interval  computation.  The  time  intervals  are  those  periods  over  which  the 
succession  of  linked  steady— state  solutions  are  obtained  for  particle  density  and  potential 
distributions.  An  analogous  parameter  is  required  for  the  algorithm  relating  ion  edge 
density  dependence,  as  well. 

For  ions,  this  parameter  can  be  calculated  readily  from  the  data  available  in  each 
code  computation;  simply  summing  the  ratio  of  spatial  grid  width  to  local  grid  average 
radial  speed,  over  the  complete  grid  range  characterizing  the  problem.  This  works 
eSectivdy,  even  though  the  ions  are  specified  with  a  square  function  energy  distribution, 
because  all  of  the  ions  go  to  or  through  the  central  core  of  the  system. 


This  is  not  the  case  for  the  electrons,  however,  as  their  radial  and  transverse  energy 
spreads,  coupled  with  the  square  distribution  used  to  model  these  spreads,  forces  them  to 
experience  (in  the  computer  code  computation)  a  wide  range  of  transit  times.  Some 
electrons  are  reflected  by  the  negative  potential  well  at  mid-range  radii,  while  others  travel 
completely  through  the  system  core  re^on.  Once  inside  the  electron-current-induced 
outer  potential  well  the  electrons  travel  inward  vrith  the  ion.  stream.  Those  at  the  top  of 
the  energy  spread  will  pass  under  the  well  (in  the  case  of  a  central  virtual  anode,  which  is 
the  potential  distribution  for  most  cases  of  interest),  and  will  exhibit  long  transit  times. 
Those  at  the  mean  energy  will  not  reach  small  radii  but  will  be  reflected  and  return  to  the 
outer  boundary  with  short  transit  times. 
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Ill  principle  it  is  possible  to  calculate  the  "average"  of  the  transit  time  for  all 
electrons  by  using  code-generated  data,  just  as  for  the  case  of  the  ions,  above.  However,  in 
practice  this  has  not  been  found  to  be  a  viable  method,  for  a  variety  of  reasons  deriving 
firom  the  numerical  mechanics  utilized  in  the  basic  computer  code  operations.  In  an  earlier 
note  the  average  dectron  transit  time  was  derived  from  simple  analytic  formulae  for 
the  transit  of  electrons  of  any  given  energy  through  potential  wdls  of  simple  shape.  These 
included  simple  power-law  potential  wells,  with  well  shape  given  by 

e0(r)  =  EJl  -  {x/Rr]  =  E^(l  -  <r>"')  (1) 

over  the  range  r^  <  r  <  R,  and  constant  at  E  »  E^  inside  the  convergence-limited  core,  at 

r  <  r  . 
o 

The  transit  time  was  also  given  for  potential  distributions  formed  of  the  simple 
power-law  well  with  a  superimposed  central  anode  potential  whose  shape  varied  linearly 
inversdy  with  radius  as 

1^0  ■" 
or 

The  transit  time  was  computed  by  integration  of  the  dectron  speed  taken  as  though 
all  electrons  were  at  the  potential  given  by  these  expressions,  from  the  integral  equation 
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where  the  integral  has  been  split  into  two  parts,  one  within  the  convergence-limited  core 
radius  (r^),  where  electron  speed  is  taken  as  constant  at  and  one  outside  this  region,  to 
the  outer  boundary  of  the  system. 

2.  INVERSE-SQUABE  CENTRAL  VIRTUAL  ANODES 

The  simple  power-law  wells  are  a  good  approximation  only  for  cases  of  very  low  ion 
and  electron  core  densities  (i.e.  for  very  low  currents).  At  high  particle  densities,  as  for 
reactor  level  conditions,  the  potential  wells  all  have  a  central  virtual  anode  formed  by  the 
convergence  of  ions  circulating  through  the  device  and  compressing  at  least  as  fast  as  the 
inverse-square  of  thdr  radial  distance  as  they  approach  the  center.  As  central  ion  density 
builds  up,  a  virtual  anode  will  be  formed  which  will,  in  turn,  slow  down  the  ions  and  lead 
to  still  greater  ion  density  rise  with  decreasing  radius.  In  the  previous  note^  the  virtual 
anode  potential  height  was  modelled  as  decreasing  inversely  with  distance  from  the  origin. 
Further  analytic  work  and  EKXL  computer  code  studies  for  relevant  cases  have  shown  that 
the  central  anode  potential  varies  more  rapidly  than  this,  proceeding  approximately  as  the 
inverse  square  of  the  distance  from  the  well  center,  thus 

1^0  -  e^.(r)l  =  (E^  -  (4a) 

or 

e^.Cr)  *  E^[l  -  (r^/r)^J  +  (4b) 

down  to  a  radius  r  ,  at  which  the  potential  well  reaches  a  minimum  value  <f>  ,  and  rises 
with  decreasing  radius  thereafter,  as  shown  in  Figure  1.  Outside  of  r^  the  electrons  slow 
down  according  to  the  variation  of  the  main  potential  (following  the  B  field).  Inside  this 
minimum-well-depth  position  the  electron  radial  speed  is  accelerated  by  the  virtual  anode 
"hill”.  The  equation  describing  the  total  potential  variation  is 
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e^t(r)  =  e^(r)-[E^-e0Jr)] 


=  EJl  -  (r/R)-]  -  EJ(r/Rrt|l  -  (e#„/E,)]  (6) 

where  e^(r)  is  the  main  system  potential  from  eq.  (1).  Differentiation  of  eq.  (5)  allows 
solution  for  the  radial  position  r^  at  which  the  net  potential  reaches  a  (negative) 
maximum  (or  well  minimum).  This  is  found  to  be  at 

(r,/il)  =  (2(t_,/R)*(EJEJ/ml‘/("‘+*)  (6a) 

or 

<r^>  =  (6b) 

where  =  (E^  —  €^^)  is  the  anode  hdght  above  the  absolute  limiting  well  bottom  E^, 
and  s=  *  1  -  (e^^/E^).  The  minimum  energy  of  electrons  at  this  point  is 

For  the  case  m  *  3,  <r^>  =  lE-3,  E^  =  2E4  eV  and  E^  =  1.1E5  eV,  for  example, 
eqs.  (6)  give  <r^>  —  0.041;  the  EKXL  code  shows  the  wdl  minimum  to  be  essentially  flat 
from  0.035  to  0.090  (mean  «  0.062).  For  a  similar  case,  vdth  the  core  at  <r  >  =:  lE-2,  it 
is  found  that  <r^>  *  0.104  and  E^  =*  309  eV  from  eq.  (7). 

The  device  can  be  operated  so  as  to  produce  any  desired  central  anode  potential 
by  adjusting  the  available  ion  injection  current  for  any  given  dectron  drive  current.  The 
EKXL  code  uses  this  central  anode  potential  as  a  target  control  parameter  in  its 
calculational  algorithm  for  adjustment  of  ion  edge  (input)  density  in  successive  time 
increments  of  computation.  Real  physical  devices  also  can  be  controlled  in  this  fashion. 
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A  closed  analytic  form  giving  an  approximate  transit  time  t^^^  for  electrons  moving 
only  radially,  with  maximum  injection  speed  v^,  across  the  system  can  be  estimated  by 
integration  over  three  terms:  one  for  the  core  region  r  <  r^,  where  the  electron  energy  is 
assumed  constant,  £  s  n£  E  /2,  at  speed  v  s  (2£  /m  one  for  r  <  r  <  r  , 
where  v(r)  =  and  one  for  the  region  r^  <  r  <  R,  where 

v(r)  =  [(2E^  +  E^^)(r/R)“/mJ®'®.  Note  that  this  integration  will  yield  an  expression 
valid  ONLY  for  maximum  energy  electrons  with  no  transverse  momentum,  and  strictly 
valid  only  for  electron  motion  under  the  full  wiffle-ball  (WB)  mode  of  operation.  Systems 
operating  at  conditions  below  the  full  WB  state  will  exhibit  slightly  longer  transit  times. 

« 

Carrying  out  these  integrations  yidds  a  general  expression  for  the  ratio 
(t.^  =  2R/v^,  the  "free-run"  transit  time)  in  terms  of  the  power  law  exponent  (m),  the 
electron  energy  £^  at  injection  and  the  (small)  energy  height  £^^  of  the  well  at  the 
maximum-well-depth  or  "turning-point"  radius  ratio  <r^>»  the  core  convergence  radius 
<r^>  and  the  anode  height  parameter  {rj).  For  the  case  of  interest  when  m  s  3,  this  is 
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The  general  behavior  of  this  transit  time  ratio  is  shown  in  Figure  2,  for  several 
values  of  the  anode  height  parameter  {rf}  and  for  several  convergence  ratios  <r^>-  Note 
that  this  algorithm  suggests  that  the  transit  time  becomes  long  compared  to  the  "free-run" 
transit  time  t,  at  small  values  of  r  and  n.  This  result  is  correct  only  for  those  electrons 
that  move  radially  with  maximum  energy,  at  the  top  of  the  injection  energy  distribution. 

It  is  clearly  non-physical  in  the  real  machine,  where  the  electrons  have  a  significant  spread 
in  energy  in  both  radial  and  transverse  directions. 
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3.  AVERAGE  TRANSIT  TIME 

To  account  for  this  it  is  necessary  to  consider  the  radial  motion  of  electrons  all 
across  their  radial  and  transverse  energy  distribution.  For  the  square  distribution  used  in 
the  EKXL  code,  an  approximate  expression  for  the  average  radial  speed  v  across  this 

avg 

distribution  is  found  to  be  given  by 

=  2/[l+(dE^.^/E,)]  (9) 

This  factor  can  account  for  the  variation  of  transit  time  of  all  those  electrons 
moving  only  radially,  over  the  function  that  describes  their  radial  energy  distribution. 

' '  However,  the  electrons  also  possess  transverse  energy  whose  "angular”  momentum  is 

conserved  in  their  recirculating  flow.  Thus  most  electrons  do  not  follow  purely  radial  paths 
but,  rather,  move  past  the  core  with  relativdy  high  speed  far  above  the  well  bottom,  at 
distances  considerably  away  &om  the  system  center. 

Exact  analysis  of  this  system  can  not  be  done  by  analytical  methods,  and  is  quite 
beyond  the  capabilities  of  any  1-dimensional  computer  code  such  as  EKXL.  An 
approximate  analysis  has  been  done,  using  simple  expansions  of  the  velocity /potential 
functions  for  the  integration  process  to  find  the  average  velocity  of  the  average  electron. 

The  correction  factor  thus  obtained  for  use  in  determing  the  average  electron  transit  time 
is  comparable  to  that  in  eq.  (9),  above,  and  is  given  by 

(10) 

With  this  the  average  transit  time  of  the  average  electron  can  be  estimated 
approximately  by  use  of  eqs.  (8, 10)  according  to 
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(t  /t.  )  f  (r  /t.  ) 

'  *vg'  in'  tm'  rtrn'  in' 


(11) 

This  is  shown  in  Figure  3  for  a  range  of  anode  height  parameter,  17,  over  a  range  of 
core  convergence  radius,  r^.  Note  that  the  average  transit  time  is  about  twice  the 
"tree-run"  transit  time  over  a  wide  range  of  parameters  of  interest  for  these  systems.  Thus 
although  the  existence  of  the  potential  well  "slows  down"  the  transit  of  electrons  passing 
through  it,  the  degree  of  this  "slowing-down"  is  not  great.  This  factor  (eq.  11)  is  used  in 
the  EKXL  code  in  the  determination  of  electron  transit  time,  for  internal  computations  of 
electron  lifetime  and  system  population. 
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Figures. 


General  potential  distribution  in  well  with  central  virtual  anode, 
=  (E^  —  e*^),  with  potential  maximum  (well  minimum)  at 
turning-point  radius,  r^ 


Transit  time  of  most  energetic  radially— moving  electrons  across  potential 
wdl  with  central  virtual  anode,  as  a  function  of  anode  height  parameter, 
V  =  and  core  convergence-limited  radius,  <r^>  =  r^/R- 


Approximate  average  transit  time  of  electrons,  averaged  over  their  initial 
radial  and  transverse  velocity  distributions,  through  potential  well  with 
central  virtual  anode,  as  a  function  of  n  and  <r  >. 
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Approximate  average  transit  time  of  electrons,  averaged  over  their  initial 
radial  and  transverse  velocity  distributions,  through  potential  well  with 
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